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A FAST FREQUENCY TRANSFORMATION TECHNIQUE 
FOR TRANSFORM AUDIO CODERS 



Technical Field 

5 

This invention is applicable in the field of multi-channel audio coders which use modified 
discrete cosine transform as a step in the compression of audio signals. 

Background An 

10 

In order to more efficiently broadcast or record audio signals, the amount of informauon 
required to represent the audio signals may be reduced. In the case of digital audio 
signals, the amount of digital informauon needed to accurately reproduce the original 
pulse code modulation (PCM) samples may be reduced by applying a digital compression 
15 algorithm, resulting in a digitally compressed representation of the original signal. The 
goal of the digital compression algorithm is to produce a digital representation of an audio 
signal which, when decoded and reproduced, sounds the same as the original signal, while 
using a minimum of digital information for the compressed or encoded representation. 

20 Recent advances in audio coding technology have led to high compression ratios while 
keeping audible degradation in the compressed signal to a minimum. These coders are 
intended for a variety of applications, including 5. 1 channel film soundo-acks. HDTV, 
laser discs and multimedia. Description of one applicable method can be found in die 
Advanced Television Systems Committee (ATSC) Standard document entitled "Digital 

25 Audio Compression (AC-3) Standard", Document A/52, 20 December, 1995. 

In the basic approach, at the encoder the time domain audio signal is first convened to the 
frequency domain using a bank of filters. The frequency domain coefficients, thus 
generated, are convened to fixed point representation. In fixed point syntax, each 
30 coefficient is represented as a mantissa and an exponent. The bulk of the compressed 
bitstream transmitted to die decoder comprises these exponents and mantissas. 
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The exponents are usually transmitted in their original form. However, each mantissa 
must be truncated to a fixed or variable number of decimal places. The number of bits to 
be used for coding each mantissa is obtained from a bit allocation algorithm which may be 
based on the masking property of the human auditory system. Lower numbers of bits 
5 result in higher compression ratios because less space is required to transmit the 
coefficients. However, diis may cause high quantization errors, leading to audible 
distonion. A good distribution of available bits to each mantissa forms the core of the 
advanced audio coders. 

10 The frequency transformation phase has one of the greatest computation requirements in a 
transform coder. Therefore, an efficient implementation of this phase can decrease the 
computation requirement of the system significantly and make real time operation of the 
encoder more easily attainable. 

15 In some encoders such as those specified in the AC-3 standard, the frequency domain 
transformation of signals is performed by the modified discrete cosine transform (MDCT). 
If directly implemented, the MDCT requires 0(N^) additions and multiplications. 
However it has been found possible to reduce the number of required operations 
significantly if the MDCT equation is able to be computed in a from that is amenable to 

20 the use of the well known Fast Fourier Transform (FFT) method of J.W. Cooley and J.W. 
Tukey (1960). A known application of an FFT to an MDCT is disclosed in, for example, 
EP-A-0564089. 

The present invention seeks to provide an alternative computation method using a Fast 
25 Fourier Transform. Moreover, the invention seeks to use a single FFT for two channels 
to achieve greater reduction in computational requirements of the system. 

Summarv of the Invention 

30 In accordance with the present invention there is provided a method for coding audio data 
comprising a sequence of digital audio samples, including the steps of: 

i) multiplying the input samples with a first trigonometric function factor to 
generate an intermediate sample sequence; 
35 ii) computing a fast Fourier transform of the intermediate sample sequence to 

generate a Fourier transform coefficient sequence; 
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iii) for each transform coefficient in the sequence, multiplying the real and 
imaginary components of the transform coefficient by respective second trigonometric 
function factors, adding the multiplied real and imaginary transform coefficient 
components to generate an addition stream coefficient, and subtracting the multiplied real 

5 and imaginary transform coefficient components to generate a subtraction stream 
coefficient; 

iv) multiplying the addition and subtraction stream coefficients with respective 
third trigonometric function factors; and 

v) subtracting the corresponding multiplied addition and subtraction stream 
10 coefficients to generate audio coded frequency domain coefficients. 

The present invention also provides a method for coding audio data, including the steps 
of: 

combining first and second sequences of digital audio samples from first and 
15 second audio channels into a single complex sample sequence; 

determining a Fourier transform coefficient sequence as defined above; 

generating first and second transform coefficient sequences by combining and/or 
differencing first and second selected transform coefficients from said Fourier transform 
coefficient sequence; and 
20 for each of the first and second transform coefficient sequences, generating audio 

coded frequency domain coefficients as defined above, so as to generate respective 
sequences of said audio coded frequency domain coefficients for the first and second audio 
channels. 

25 The present invention also provides a method for coding audio data including the steps of: 

obtaining at least one input sequence of digital audio samples; 

pre-processing the input sequence samples including applying a pre-multiplication 
factor to obtain modified input sequence samples; 

transforming the modified input sequence samples into a transform coefficient 
30 sequence utilising a fast Fourier transform; and 

post-processing the sequence of transform coefficients including applying first post 
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muliiplicaiion factors to the real and imaginary coefficient components, differencing and 
combining the post-multiplied real and imaginary components, applying second post- 
multiplication factors to the difference and combination results, and differencing to obtain 
a sequence of modified discrete cosine transform coefficients representing said input 
5 sequence of digital audio samples. 

The present invention also provides a method for coding audio data including the steps of: 

obtaining first and second input sequences of digital audio samples corresponding 
to respective first and second audio channels; 

10 combining the first and second input sequences of digital audio samples into a 

single complex input sample sequence; 

pre-processing the complex input sequence samples including applying a pre- 
multiplication factor to obtain modified complex ir^ut sequence samples; 

transforming the modified complex input sequence samples into a complex 

15 transform coefficient sequence utilising a fast Fourier transform; and 

post-processing the sequence of complex transform coefficients to obtain first and 
second sequences of audio coded frequency domain coefficients corresponding to die first 
and second audio channels including, for each corresponding frequency domain coefficient 
in the first and second sequences, selecting first and second complex transform coefficients 

20 from said sequence of complex o-ansform coefficients, combining the first complex 
transform coefficient and the complex conjugate of the second complex transform 
coefficient for said first channel and differencing the first complex transform coefficient 
and the complex conjugate of the second complex transform coefficient for said second 
charmel, and applying respective post-multiplication factors to the combination and 

25 difference to obtain said audio coded frequency domain coefficients corresponding to the 
first and second audio channels. 

The present invention further provides A method for coding audio data including the steps 
of: 

30 obtaining first and second input sequences of digital audio samples x[n], y[n] 

corresponding to respective first and second audio channels; 
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combining the first and second input sequences of digital audio samples into a 
single complex input sample sequence z[/z], where z[n] = x[n] -h Jy[n] ; 

pre-processing the complex input sequence samples including applying a pre- 
multiplication factor cos{m/N) -h Jsin(m/N) to obtain modified complex input sequence 
5 samples, where N is the number of audio samples in each of the first and second input 
sequences and n = 0,....,(A^-1); 

transforming the modified complex input sequence samples into a complex 
transform coefficient sequence utilising a fast Fourier transform, wherein k = 
0 ,(M2-l);and 

0 post-processing the sequence of complex n^form coefficients to obtain first and 

second sequences of audio coded frequency domain coefficients corresponding to the first 
and second audio channels Xj,, according to: 



= cosy * (^^cos(7t(A-^l/2)/A0-^j^.sin(7t(>t+l/2)///) 

- siny * (^^^sin(n(Ar+I/2)/A^+gj^.cos(7t(;t+I/2)/A0 

= cosy * (g'jfc^,cos(7t(A^I/2)/A0-g'jtwS"^(^(*^^/2)/A0 

- siny * C?'jfc.^sin(7t(/r+l/2)/A0*g'^.cos(7t(A:+l/2)/A0 

Gt is a transfonn coefficient sequence for die first channel; 
G'jt is a transform coefficient sequence for the second channel; 
gi.^r and g^- are the real and imaginary transform coefficient components of G^; 
g and g\- are the real and imaginary transform coefficient components of 
T ^.jt.^ is the complex conjugate of Z^sr-t-/; and 
yW = 7r(2A+l)/4. 

modified discrete cosine transform equation can be expressed as 



15 where 



20 
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^1- = H * cos(27i*(2/7+l)*(2A-^l)/4iV + 7t *(2A:+l)/4) ^=0...(M2-1) 

where xfnj is the input sequence for a channel and N is the transform length. 

Instead of evaluating A; in the form given above it could be computed as 

X^, = cosY*(^^cos(7t(fc+l/2)/A0-5^^,sin(7r(;t+l/2yA0) 
-siny ^(gf^^sm(Tz(k^ l/2yN)^g^^cos(n(k^ \/2)/N)) 
Sk^^kj ^ S*(-y^^ of real numbers) 



where G, = g^^ ^ jg^^ = . The symbol j represents the 

n=0 



5 imaginary number . The expression Yl {x[n\e^"'^)*e^'^'''^^ is obtained from 



the well known FFT method, by first using transformation a:* [/2]=x[rz] * e^"^^ and then 
computing the FFT - ^ t^] ^e^'^"^^ 



For a two channel approach, a complex variable z[n] = xM^e^'^^ 4- jy[n]*e^""^ is 
10 defined, where x[n] and y[n] are sample sequence for the two channels and e^'™^ 

represents the pre-multiplication factor. Using FFT approach, the frequency coefficient 
for the variable z[n] is computed. From Z,, the value = (Zj, + Zr^.,,J/2 and G; ={Zj, - 
2"//.jk.;)/2/\ required to compute the final MDCT for each channel, respectively, is 
calculated. 

15 

If either or both the channels require short length u-ansformers, two shon transforms are 
taken using the above approach. If neither need shon transform, a single long transform 
is used. As an additional step in reducing computation, the windowing function can be 
combined with the pre-processing stage. 



r 
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Brief Description of the Drawings 

The invention is described in detail hereinafter, by way of example only, with reference to 
preferred embodiments thereof and with aid of the accompanying drawings, wherein: 
5 Figure 1 is a diagrammatic representation of a stream of audio data and the 

substructure arrangement thereof; 

Figure 2 is a functional block diagram of a digital audio encoder; 

Figure 3 is a functional block diagram of a system for encoding a single audio 
charmel; and 

10 Figure 4 is a functional block diagram of a system for encoding a pair of audio 

channels. 

Detailed Description of the Preferred Embodiments 

15 The above mentioned Advanced Television Systems Committee (ATSC) Standard 

document entitled "Digital Audio Compression (AC-3) Standard" (Document A/52, 20 
December, 1995) describes methods for encoding and decoding audio signals, and is 
hereby expressly incorporated herein by reference. 

20 In general, the input to an audio coder comprises a stream of digitised samples of the time 
domain analog signal. For a multi-channel encoder the stream consists of interleaved 
samples for each channel. The input stream is sectioned into blocks, each block 
containing N consecutive samples of each channel (see Fig. 1). Thus within a block the N 
samples of a channel fomi a sequence {;c[0], x[l], x[2], ^A'"-!]}. 

25 

The time domain samples are next convened to the frequency domain using an analysis 
filter bank (see Fig. 2). The frequency domain coefficients, thus generated, form a 
coefficient set which can be identified as (Xo. X,, X,, X ^.r^^i)- Since the signal is real 
only the first N/2 frequency components are considered. Here X^ is the lowest frequency 
30 (DC) component while X^^^,, is the highest frequency component of die signal. 
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Audio compression essentially entails finding how much of the information in the set (Aq, 
) is necessary to reproduce the original analog signal at the decoder with 
minimal audible distortion, 

5 The coefficient set is normally converted into floating point format, where each coefficient 
is represented by an exponent and mantissa. The exponent set is usually transmitted in its 
original tbnm. However, die mantissa is truncated to a fixed or variable number of 
decimal places. The value of number of bits for coding a mantissa is usually obtained 
from a bit allocation algorithm which for advanced psychoacoustic coders may be based 
10 on the masking property of the himian auditory system. A low number of bits results in 
high compression ratio because less space is required to transmit die coefficients. 
However this causes very high quantization error leading to audible distonion. A good 
distribution of available bits to each mantissa forms die core of the most advanced 
encoders. 

15 

In some encoders such as the AC-3, the frequency domain transformation of signals is 
performed by the (MDCT) modified discrete cosine transform (Eq. 1). 

= * cos(2k*(2/7-*-1)*(2/:+1)/4A^ * t:*(2*+1)/4) >t=0...(A72-l) Eq 1 

If direcdy implemented in the form given above, the MDCT requires O(l^) additions and 
multiplications. 

20 

Sinele Channel FFT 

It is possible to reduce the number of required operations significantly if one is able to 
25 evaluate Eq. 1 using the well known Fast Fourier Transform method of J.W. Cooley and 
J. W. Tukey (1960). The general Discrete Fourier Transform (DFT) is given below (Eq. 
2). It requires O(N^) complex additions and multiplications. By using the Fast Fourier 
Transform method the DFT in Eq. 2 can be computed with 0{Nlog2N) operations only. 
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^1-= E M"]*^^''"*^^) A:=0..^-1 Eq. 2 

rt=0 

Herey is the symbol for imaginary number, i.e. y = yAL 

Although it may not be immediately apparent how Eq. 1 can be transformed to Eq. 2, a 
careful analysis shows that this is indeed possible. To simplify Eq. 1, two functions can 
5 be defined 

cL{n,k) = 271(2/2 + l)(2Jt+l)/4N - Eq. 3 

Y(*) = ^(2/:+l)/4 Eq. 4 

Then, using these functions, Eq. 1 can be rewritten as 



^A: = H 4^]^cos(a(/2,A)+YW) Ea 5 

= ^M*(cosa(/2,A)cosY(Ar)-sina(/i,A:)sinY(>t)) Eq. 6 



/r=0 



10 In Eq. 6 the trigonometric equality, cos(a+f7) = cosa cos6-sina sinZ? is used for 

simplification. Furthermore, since the function Y(k) is not dependant on variable it can 
be brought outside the summation expression to give 

= cosyik) J2 xln]*cosa{n,k) - sinYW Z x[n]^sina(n,k) 

= TjCOSyC*) - T^sinyik) ^ 

Where = E x[n]*cosa(n,k) and 7^ = £ *sina(n,Jt) 



The two terms, and T^. can now be evaluated separately. Using Euler's identity e*** = 
15 cose -hysinO, we can express: 

cosa(/i,/:)=(e'**'^*'* + e'^"*"'*'')/2 
and sina(/z,/:) ={e'''*''-''*-e'j"*"**'*V2y . 
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Therefore we can rewrite the term Tj as 

r, = ^ x[/7].(e>«+e-^«)/2 = l/2( x[n]*e^ + J2 ^W-^"''") 

= 1/2 - A^) 



where = 4^^]*^^" and = 52 4^]*^"'^'" 

n=0 * rt=0 



Similarly 



^2 = E xM*(e^''-e-^'*)/2 = l/2/( ^ x[n]*e>« - E^W*^'''') 



The term can thus be evaluated from Eq. 8 and Eq. 9 

ff=0 



Eq. 10 (T 



5 If a complex variable is defined as: 

x't/2]=x[n]*e*'^^ Eq. 11 

then Eq. 10 is simply: 



Eq. 12 



where G^. = x^[w] ^e^^''"^^ 

n=0 
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The complex term = gk,r+gkj» where g^,^ and ^ e 3t (set of real numbers) in Eq. 12 is 
essentially the same as in Eq. 2. Therefore the FFT approach can be used to evaluate 
Gi,. This brings down computation from O(N^) to O(NlogN). Similarly, the second term 
A2 in Eq. 8 and Eq. 9 can be evaluated 



Eq. 1: 



5 where = Yl {x[n]^e'^"^^)*e'^^'''^^ 



Note diat G/ is actually the conplex conjugate of G^ which was obtained by Eq. 12. That 
is, if = g,.,+gfc,i, where g,,, and g^ e 31 as defined earlier, then G/ = g^^, - yg^ ^. 
Therefore G^' in Eq. 13 does not need to be computed again, and the result from Eq. 12 
can be re-used. That is, only one FFT needs to be confuted for the evaluation of T^, 
10 The result of Eq. 8 to Eq. 13 is thus 



= l/2(e^(^^^^>^ G^^e 'M^^^^W c?;) Eq. 14 



Next, the term r, can be analysed 



7; = E x[nMe^''-e'n/2j 



Finally, after simplifications of Eq. 7, 14 and 15 
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= cosY(^) l/2(e^<^*'^>^ G^^e-^^^^^"^^^ G^) 

= cosy * (gj^cos(n(k^l/2)/N)'gj^^sm(Tz(k-^\/2)/M) 
- siny * (g^^sin(Tz(k^l/2)/N)^g^cos('K{k-i'l/2)/N) 
= cosy - siny * 



The term G^ = g^^ + y^^. is computed in O(NlogN) operation by use of FFT algorithms. 

The additional operation outlined in Eq. 16 to extract the final Xj, is only of order 0(N). 

Therefore the MDCT can now be computed in 0(N\QgJSI) time. The operations required ^ 

to obtain the MDCT are illustrated in Fig. 3. 

5 

Combining Two Channels into Single FFT 

Suppose the multi-channel encoder is required to process m audio channels. Instead of 
computing an FFT for each channel as described in the previous section, it is possible to 
10 further reduce the computational requirement of the coder by combining two channels and 
using a single FFT only. In effect, instead of m FFTs only m/2 FFTS need to be 
computed. 

If the input sequence are real numbers then it is known that DFT for any two channels can 
15 be computed with only one FFT block by considering the input as a complex number. 
The real pan is formed from the sequence for any one channel and the unaginary part is 
from data of another channel. After the Fourier Transform is computed for the resulting 
complex variable, the resulting transform for each channel can be easily retrieved. 

20 However, in the present case the input data to the FFT block is actually a complex number 
(formed by multiplying the real data by complex variable e^"^'^. In this case, there is no 
straightforward way of retrieving the frequency transform after having combined two 
channels. However, using some processing after the FFT one can still compute the DFT 
of two channel using a single FFT block. 



25 
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Let {j:[0]^[1]^[2],...^[N-1]} be input samples of the first channel and 
{yL0],y[l].y[2],, . .yU^-l]} be the samples for the second channel. As described above, the 

frequency coefficients = J2 ^W^^'"^ * ^j2iink/f/ ^2 and 13) must be 

rt=0 

obtained for the first channel; and similarly, for the second channel 



n-N'\ 



nnk/N 



rt=0 



Defming complex variable z[«] = jcC/zjV"^^ + y>[n]*e'"°'^ Eq. 17 

and computing its DFT using the FFT method, yields 

"^0 Eq. 18 

Now substituting N-k for k in the above expression. 



Eq. 19 



Since e'-"*' - 1, /i e I (the set of integers), the term d^""" vanishes in the above expression. 
10 Taking the complex conjugate of Z^^-t- 

= "e' (:cM-y>M).e-^^-(-^^^'^^VN 
= IT (x[A7]-y>[/,])*e^-''"<*-"2>'V 

Eq. 20 
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Using Eq. 18 and 20, separate expressions for and G\ are required. In a simple case 
the conjugates in Eq. 18 and 20 should add and subtract to give the required expressions. 
However in this instance that is not the case. But, substituting A''-* by N-k-I in Eq. 18, the 
following is obtained 

= "f: ixin]-jy[n])^e^-'-"<'^'^y^ Eq. 21 

5 Now the term e'-"'^*^"*'^^'^ is conunon in both Eq. 17 and 19, and it is possible to isolate. 



ff*0 

= 2G, 



Similarly, 



That is 



= (Zj^+Z;.jt.,)/2 Jt=0..JW2-l Eq. 22 



and 



10 



k=0..M2-\ 



Eq. 23 
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From the expression from Eq. 22 and 23 into Eq. 16, the MDCT for each channel is 
obtained. The overaJl process is illustrated in Fig. 4. 

Transform Length Adjustment Technique 

5 

The frequency transform length N is decided by the encoder based on temporal and 
spectraJ resolution requirements. The input signal is usually analysed with a high 
frequency bandpass filter to detect the presence of transients. This information is used to 
adjust the block length, restricting quantization noise associated with the transient within a 
0 small temporal region about die transient, avoiding temporal masking. Thus, if transient 
is detected in a channel, two short transform of length N/2 each are taken. In the absence 
of transient, a single long nransform of length N is used, thus providing higher spectral 
resolution. 



From the method described in the previous secdon for computing MDCT for two channels 
using a single FFT block, it is evident that the transform length for the two paired 
channels must be the same. Therefore, pairing for the transformation phase much be such 
that channels with identical u-ansform length are grouped together. 

It is however possible that not all channels can be paired with such convenience. Assume 
that the total number of channels are an even number (if not, take a single FFT for one 
channel and the rest form an even group). Suppose out of the m channels, / need long 
transform and therefore m-l require shon transform. 



25 It / is an even number, then since the total is even, it follows that l-m is also even. In this 
case, from the / channels that need long transform, 1/2 pairs are formed and for each of 
the 1/2 pairs a single FFT is computed to estimate the MDCT for the original paired 
channels. Similarly, the l-m channels are paired to form (l-m)/2 pairs and for the (l-m)/2 
pairs two short FFTs are computed. 

30 



consider the case when W = 2r + 1 is an odd number. Therefore m - / = 2s + 1 is 



m 
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also an odd number. The 2r channels requiring long transform are paired together to form 
r pairs and then 2r transforms are computed using r FFTs only. Similarly, for the 2s 
channels s pairs are formed. What remains is one channel requiring long transform and 
another requiring two short tfansforms. Both of diese channels are paired togedier and 
5 two short FFTs are computed to derive the MDCT. 



The rationale for constraining die long transform to two short ones is as follows. A short 
transform is required for restricting quantization noise associated with the transient within 
a small temporal region about the transient, avoiding temporal masking. A long transform 
10 gives slight bener frequency resolution but the error is not much compared to the case ^ 
when in the presence of transient a long transform is utilised. Forcing a long transform 
onto a channel in the presence of transient leads to greater distortion in the final produced 
music. This conjecture was proven true by experimental studies on benchmark music 
streams. 

15 

Combining Windowing with pre-processing 

Before the time domain signal x[n] is transformed to the frequency domain, a windowing 
ftmction is usually applied. Thus, if the sampled signal is p[n] then the sequence that is 
20 applied to the frequency transformation block isx[n] = p[n] * w[n], where w[n] is the 
windowing function. From die previous sections we noted that before die FFT is 
computed for a block a pre-processing is performed as given in Eq. 11 (reproduced below 
for convenience). Thus 

x'[n] = x[n] * e'""'^ 

= (pW* wCrt])* (cosTmm + J sin rm/N) 

= * ((w[«] * CQSJOi/N) -h j{w[n]*sm7m/N)) Eq. 24 

From Eq. 24 we note that the windowing function can be combined wiUi the cosine and 
30 sine multiplication required in Eq. 1 1. This brings down the computation even funher 
since the sine and cosine are usually implemented in a real time system as table-lookup. If 
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two tables are constructed as defined below 

rcos[n] = w[n] * cos(7m/iV) 
rsin[72] = w[n]* sin(;r7z/A0 

5 

then Eq. 1 1 can be rewrinen as 

x'[n] = (j)[n] * rcosM) + J(p[n]* r sin [n]) Eq. 25 

0 Although the invention has been described herein primarily in terms of its madiematical 
derivation and application, and the procedures required for implementation, it will be 
readily recognised by those skilled in the art diat the procedures described can be 
implemented by means of any desired computational apparatus. For example, the 
invention may be embodied in computer software operating on general purpose computing 

5 equipment, or may be embodied in purpose built circuitry or contained in microcode or 
the like in an integrated circuit or set of integrated circuits. 



20 



The foregoing detailed description of embodiments of die invention has been presented by 
way of example only, and is not intended to be considered limiting to the invention as 
defined in the claims appended hereto. 
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Glossarv of Equations: 
MDCT 

= Yl 4"] • cos(27r*(2n+l)*(2Xr+l)/4// + 7t •(2jt+l)/4) A:=0...(AV2-1) 

(T=0 

= COSY * (Sk^cos(Tz(k+l/2)/N)-g,^sin(n(k->-l/2yN) 
- siny • (s,^sm(n(k*l/2)/N)*g^cos(Ti(k*V2yN) 
= r,cosY(-t) - T^smy(k) 



^1 = IT x[n]*cosain,k) = *sina(/7,Ar) 

= l/2(^, * A^) = l/2J(A^ - A^) 



"=0 11=0 



G^. = (Zt*Z;.^.,)/2 k=0..Nl2-\ 
10 = (2t-Z;-t.,y2/ k=0..M2-\ 



a{n,k) = 2i;(2n+l)(2A:+l)/4N 
YW = TC(2it+l)/4 



